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Abstract 

We study the evolution of the size - stellar mass relation for a large spectroscopic sample of galaxies 
in the GOODs North field up to z ~ 3.5. The sizes of the galaxies are measured from if s -band 
images (corresponding to rest-frame optical/NIR) from the Subaru 8m telescope. We reproduce earlier 
results based on photometric redshifts that the sizes of galaxies at a given mass evolve with redshift. 
Specifically, we compare sizes of UV-bright galaxies at a range of redshifts: Lyman break galaxies 
(LBGs) selected through the U-drop technique (z ~ 2.5 — 3.5), BM/BX galaxies at z ~ 1.5 — 2.5, and 
GALEX LBGs at low redshift [z ~ 0.6 — 1.5). The median sizes of these UV-bright galaxies evolve as 
(1 + z)- 1 11±013 between z ~ 0.5-3.5. The UV-bright galaxies are significantly larger than quiescent 
galaxies at the same mass and redshift by 0.45 ±0.09 dex. We also verify the correlation between color 
and stellar mass density of galaxies to high redshifts. The sizes of sub-mm galaxies in the same field 
are measured and compared with BM/BX galaxies. We find that median half-light radii of SMGs is 
2.90 ± 0.45 kpc and there is little difference in their size distribution to the UV-bright star forming 
galaxies. 

Subject headings: galaxies: evolution - galaxies: high redshift - galaxies: structure 



1. INTRODUCTION 



Recent studies provide evidence that sizes of galaxies 
at high redshifts were smaller in compariso n with galax- 
ies of similar mass in the l oc al universe ( e.g. Dadd i et al. | 
20051: iTruiillo et al.ll2006alfb1:lZirm et alj|2007t iToft et all 
2007t Ivan Dokkum et all I2008at iFranx et all 120081: 



Willi ams et all 120101; IToft et al.l 120091: IDamianov et all 
20091 : iCimatti et all 120081 : iLonghetti et al.l 12007ft . These 



studies pose several questions about the evolution of 
properties of galaxies and the build up of galaxy stel- 
lar mass with time. Several physical processes are pro- 
posed to explain the growth of galaxies with time, such as 
galaxy (major or minor) mergers (e.g. iKhochfar fc Silkl 
2006. 120091 : Bell et al.l l2006: Na ab et al.lboO^ gas accre- 
tion in outer regions and star formation (e.g. lDekel et ahl 
120091 : lElmegreen et al.l 120081 ). However, each of these 
mechanisms will affect differently the growth of mass 
and size of galaxies. For example, major dry mergers 
could produce descendants with larger than observed 
stellar masses for early-type galaxies. Therefore, some 
authors (e.g. Ivan Dokkum et ail [2010l: IB czan son et al.l 
l2009HNaab et al.ll2009l ; Hopkins et al.ll2009ft suggest mi- 
nor mergers and envelope growth by satellite accretion 
to be a more likely mechanism for building up quiescent 
galaxies. In general, plausible models have to reproduce 
galaxy sizes and stellar mass densities seen today and 
bring the high redshift galaxies to the local stellar mass- 
size relation. Therefore, exploring the observed evolution 
of galaxy sizes is essential to constrain these galaxy for- 
mation and evolution models. 
However, most of the size studies at high redshift were 
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based on photometric redshifts. Although techniques 
to construct galaxy spectral energy distributions (SEDs) 
from multiwavelength observations and photometric red- 
shift measurements are improved over the past few years, 
photometric redshifts still have worse uncertainties for 
star forming galaxies due to the lack of a strong 4000A 
break. The large uncertainties and possible systematic 
biases, due to the reliance on photometric redshifts can 
seriously affect the masses and star-formation rate mea- 
surements and thus the final results. Spectroscopic ob- 
servations of high-rcdshift galaxies and making complete 
magnitude limited samples require large amounts of ob- 
serving time. Nevertheless, the advent of 8- 10m class 
telescopes with multi object spectroscopy has brought 
about much larger high-z spectroscopic samples than 
were previously available. 

Spectroscopic redshifts are relatively easier to obtain 
for star forming galaxies (compare to quiescent galax- 
ies) due to emission lines. Previous results on the size 
evolution of this population (with photometric redshifts) 
are interesting since some of the usual mechanisms as- 
sumed for size growth, e.g., gas-poor ("dry") mergers do 
not apply to gas-rich starburst galaxies since substan- 
tial amounts of gas is required for these galaxies to un- 
dergo significant star formation. Therefore, studies of 
the structure of star forming galaxies with secure red- 
shifts can provide strong constraints on size evolution 
models. 

The actively star forming galaxies at high redshift con- 
sists of two main populations: relatively unobscured UV- 
bright galaxies and dusty red starbursts, the most ex- 
treme of which are detected as sub-millimetre galaxies 
(SMGs). The UV-bright galaxies are selected to have 
strong rest-frame UV. In the past few years, many sim- 
ple photometric techniques (e.g., "[/-dropout") were de- 
sign ed to select these galaxies at different redsh ift ranges 
(e.g. lSteidel et al.|[2003l: lAdelberger et al.ll200l . 
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The SMGs are among the most massive, luminous, 
and vigorously star forming galaxies at high redshift 
that are heavily obscured by dust (e.g. iHughes et al.l 
HoMlSmail et al.ll2lMIChaoman et al.l2003Ll2005D . The 
physical process driving these highly luminous galaxies 
(star- formation, AGN or combination of these two) is still 
uncertain. The evolution of SMGs and their relation to 
the local galaxies are also not yet understood: for in- 
stance, whether or not SMGs are the progenitors of local 
ellipt i cal galaxies (e.g.lBlain et af1l2004t iSwinbank et al.l 
120061 : iTacconi et al.l 120081 ) . Their position on the size- 
mass plane may therefore provide clues about their rela- 
tion to other galaxy populations. 

To verify previous results based on photometric 
redshifts, here we study the mass and size evolution of 
a large sample of UV-bright and submillimetre galaxies 
in GOODS-North with secure spectroscopic redshifts. 
The structure of this paper is as follows. In section 2, 
we review the data. In section 3 & 4, we describe size 
and mass determinations of galaxies. Our sub-samples 
selection are described in section 5. Finally, in section 6 
we present our results and investigate the size evolution 
and stellar mass-size relation for our samples. We 
summarize and discuss our results in section 7. The 
cosmological parameters adopted throughout this pa- 
per are Cl m = 0.3, il\ = 0.7 and Hq = 70 km s _1 Mpc -1 . 



2. DESCRIPTION OF DATA 

The sample of galaxies used here is based on the 
most complete spectroscopic catalog of galaxies in 
the GOODS (Great Obs ervatories Origi n Deep Sur- 
vey dGiavalisco et all 120041) ) North field by IBarger et al.l 
(2008). This catalog gives a compilation of all spec- 
troscopic observations carried out in this field (e.g . 
Cowie et al.l [200l IReddv et al.|[200l IWirth et all [2001 
Cohenl I200ll: ICohen et al.l 120001) . where each galaxy 
sample was s electe d in a different way. In addition, 



iBarger et al.l (|2008l ) performed spectroscopic observa- 
tions for certain subsamples. The catalog includes 2907 
sources including stars and is restricted to sources with 
K SlA B < 24.5 or F850LP AB < 26. There are 2362 
sources with z < 1.6 and 327 sources between z = 1.6 
and 3.5. Most redshifts for galaxies with z > 1.6 come 

from IReddv eVaTl (I2006T) . 

This catalog (jBarger et all I2008D comprises optical 
photometric data in the F435W, F606W, F775W and 
F850LP passbands tak en from HST Advanc ed Cam- 
era for Surveys (ACS) (jGiavalisco et al.ll2004D a nd th e 
U-band magnitude is taken from ICapak eTaLl (l200l . 
The near-infrared iVj-band magnitude measured from 
WIRCam images from CFHT. Details of near-IR obser- 
vations, data re duction and g enera ting catalog are de- 
scribed more in IBarger et all ([2008I ). X-ray soft (0.5 - 
2KeV) & hard (2 — 8KeV) luminosities are also pro- 
vided for many sources in the catalog. In Barge r et al.l 
(2008) sources with X-ray luminosities above 10 ergs 
s _1 in either soft or hard band defined as AGNs. For a 
fraction of sources (<~ 36%) near-ultraviolet (NUV) and 
far-ultraviolet (FUV) magnitudes from the UV imaging 
survey performed by GALEX mission were also provided. 
Note that sources indicated as stars in the catalog are ex- 
cluded from our analysis. 



3. SIZES 
3.1. Size measurements 

Size determination in the observed near-infrared (rest- 
frame optical at high z) is more robust in tracing the 
distribution of stellar mass than the observed opti- 
cal (rest-frame UV at high z) since the sizes at rest- 
frame UV can be strongly affected by dust and star- 
formation. Therefore, we use _ftT s -band images from 
Subaru MOIRCS Deep Survey (MODS) in the GOOD S 
North field (jKaiisawa et al.l l200l lOuchi et al.l l200l . 
The imaging observations are performed for J and K s 
bands over ~ 112 arcmin 2 area of the GOODS-N field 
with a pi xel scale of 0.12 arcsec . We used the images re- 
duced bv iBouwens et alj (|2008D . The resulting FWHM 
of the K s images is ~ 0.5". The deepest data cover ~ 28 
arcmin 2 reaching a 5 a of 25.4 AB mag in K s band while 
the other regions are ~ 1.3 mag shallower. More de- 
tails of the images an d data reduction can be found in 
IBouwens et"aT] (|2008f >. 

Sizes of galaxies were es timated by using the GAL- 
FI T code (iPeng et al.ir2 002) in a similar procedure used 
bv lTruiillo et al.1 (j2007ft and iWiiliams et alj pOlCft . For 
each galaxy, a square postage stamp of 120 pixels (14.4") 
around each galaxy was made and we use a mask to ex- 
clude neighbouring galaxies from the fit. A range of 
Sersic (1968) profile models convolved with the point- 
spread function (PSF) of the image are fitted to the K s 
images of galaxies. The convolved models for each object 
were compared with the galaxy surface brightness dis- 
tribution and the best fit model determined using min- 
imized x 2 of the fit. The PSF used by GALFIT was 
taken from the median of the unsaturated stars over the 
entire MOIRCS K s images. We perform a test to check 
the PSF of the image by m easuring sizes of a n umber of 
stars in the field listed in the Barger et al. (2008) catalog. 
The derived sizes of all the stars are found to be <C 0.02 
pixel (effectively zero). Therefore, the PSF is a good ap- 
proximation of a point source. We note that within our 
K s band images, there is a shallower region with seeing 
of 0.68". We followed a similar procedure with an ap- 
propriate PSF (which is made from the median of stars 
in this region) to measure sizes of galaxies in this region 
separately. We also verified that the derived sizes of the 
shallower region are reliable by its appropriate PSF. 

We determine the circularized effective radius r e = 
a e y(l — e) , from the half-light radius along the semima- 
jor axis a e and ellipticity (e) as output by GALFIT. This 
removes the effects of ellipticity. We allowed the Sersic 
index n (which measures the shape of the surface bright- 
ness profile of galaxy) to vary between 0.5 and 5 and 
the effective radius between 0.01 and 60 pixels (0.0012" 
and 7.2"). Initial guesses for the effective radius r e , el- 
lipticity and position angle were taken from the SExtrac- 
tor catalog and ma gnitude was taken from the original 
(Barg er" et al.ll2008h catalog, and we set the Sersic index 
to 2 initially. 

Our results show that the median Sersic index mea- 
sured for all galaxies is 2.0 where 50% of measurements 
lies between 4.0 and 1.0. For galaxies with stellar masses 
between 10 10 -10 11 Mq the median Sersic index is 2.4 
where 50% lies between 4.6 and 1.2. 

Due to possibility of color gradients, it is best to use the 
same rest-frame band for measuring sizes of galaxies at 
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Fig. 1. — Left panel: The points are the median of the relative difference between the recovered and input sizes versus magnitude based 
on our simulations. Right panel: The same comparison but for the relative difference between output and input Sersic index. The error 
bars illustrate the 68% scatter. The random uncertainties of recovered sizes and Sersic index increase with magnitude. 



all redshifts. However, only deep K imaging is available, 
so there is a possib ility of systematic effects with redshift. 
iFranx et al.l ((2008) show that such systematic effects are 
small and will not significantly affect the results. As 
an approximation, we corrected sizes to the rest-frame g 
band u sing star forming galaxies from studies of CDF- 
South (jFranx et al.l 12008ft . We derive the best linear fit 
to the median of the ratio of sizes in K band and rest- 
frame g band as a function of redshift and apply this 
fit (r et k-bcmd/r e , g -band = 0.15z + 0.64) to galaxies with 
z < 2 in our sample. 

3.2. Simulations 

To gauge the accuracy of the measured sizes and re- 
liability of our results we perform a realistic simula- 
tion. About 8500 Sersic profiles were generated with 
uniformly-distributed random parameters in the ranges 
of 19 < K AB < 25, 0.1" < r e < 2", 0.5 < n < 5 and 
< e < 0.8. The mock galaxies were then convolved with 
the PSF of the image. Finally we added these galaxies 
to 14.4" blank-sky postage stamps randomly taken from 
our K s band image. The structural parameters of the 
model galaxies were then measured in a manner identi- 
cal to that used for the actual images. 

The results of these simulations are shown in Figure 1 
and 2. Left panel of Figure [1] shows that the scatter in 
the recovered sizes increase with magnitude. The ran- 
dom uncertainties increase significantly after Kab = 23. 
However, the systematic errors are very small (< 5%) 
even at the faintest magnitudes. Moreover, the simu- 
lation shows that the systematic errors on recovery of 
Sersic index are also very small and random uncertain- 
ties increase above K =23 (right panel of Figure[IJ. From 
this result, we limit the sample studied in this paper to 
galaxies with K < 23 mag. We note that this limit is 
for the deeper region in our K band images. Running 
a separate simulation for the shallower region, we find 
that the derived sizes of galaxies with K < 22.7 mag are 
reliable. Hence, we restricted the sample in this region 
to this slightly brighter magnitude limit. 



We have also explored how the recovery of size depends 
on the size itself. As shown in Figure [2] the uncertainties 
in retrieving sizes depend on both magnitude and size. 
Systematic offsets at all magnitudes appear to be neg- 
ligible except for r e > 1.8 arcsec (which corresponds to 
14.4 kpc at z ~ 1) at the faintest magnitudes; however, 
galaxies this large and faint have not been seen at any 
redshift. 

The random uncertainties increase with size in all mag- 
nitude bins. For objects with Kab < 20 mag (top panel) 
the random uncertainties in size recovery are < 5%. 
However, for galaxies with 21 < Kab < 22 the increas- 
ing of random uncertainties with sizes is significant. The 
increase in random uncertainties at larger sizes could be 
due to decreasing of surface brightness. However, since 
we are mainly concerned about the overall sample prop- 
erties, the lack of systematic errors is more important. 

4. STELLAR MASS ESTIMATES 

The stellar masses of galaxies were measured with the 
Fittin g and Assessment of Synthetic Templ ates(FAST) 
code (iKriek et all [2009bh . All fluxes from IBarger et all 
(2008) including four optical bands from ACS, the K s 
infrared and U band were used to find the best-fit galaxy 
template SED to the broad band photometry using a 
X 2 minimization procedure. iBruzual fc Chariot! (2003) 
stellar population evolution models with exponentially 
declining star formation histories (with r ranging from 
10 7 ~ 10 yr) wer e used to fit the SEDs. We use the 
ISafpeterl (|1955l ) initial mass function (IMF) and solar 
metallicity and the extinction Ay was allowed to vary 
between and 3. Redshift of galaxies were fixed to 
their spectroscopic redshift provided by IBarger et all 
(2008) catalog. Masses were then shifted by —0.2 dex 
for consistency wi th the z ~ SP SS masses (which were 
calculated using a lKroupal (|2001ft IMF). 



5. SUB-SAMPLES 
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Fig. 2. — Relative difference between the input and measured 
sizes of galaxies in our simulations as a function of size for mock 
galaxies for galaxies with K^g < 20 mag (top panel), 20 < Kab < 
21 (middle panel) and 21 < Kab < 22 (lower panel). The rela- 
tive errors between the input and measured sizes of galaxies as a 
function of sizes depends on the magnitude. 

The ability to identify and study distant galaxies has 
improved dramatically during the last two decades. Var- 
ious selection criteria have been designed to select high 
rcdshift galaxies through their observed colors. In this 
paper we use such criteria to select samples of UV-bright 
galaxies in redshift ranges of 0.6 < z < 1.4, 1.4 < z < 2.5 
and 2.7 < z < 3.5. 

Our sam ple of Lyman break galaxies (2.7 < z < 3.5) is 
taken from lReddv et all (|2006l ). Their LBGs candidates 
were originally preselected by the "C" , "D" and "MD" 
criteria whi ch use regions of ( U n — G) versus (G — R) 
color space (jSteidel et al.ll2003l ). Limiting our sample to 
Kab < 23, we were only able to include 5 LBGs in our 
analysis. 

At z ~ 2, B M/BX galaxies also are taken 
from iReddv et al l (120061). The BM/B X criteria 
(|Adelberger et all 120041: ISteidel et al.1 [2004) were de- 



Sample 


Redshift 


log(Ah/M e ) 


No. of Sources 


GALEX/LBG 


0.6 < z < 1.4 


8.8 - 11.0 


105 


BM/BX 


1.4 < 2 < 2.7 


9.8 - 10.8 


11 


LBG 


2.7 < z < 3.5 


10.2 - 10.8 


5 


SMG 


0.5 < z < 3. 


10.0 - 11.7 


11 


sBzK 


1.4 < z < 3.2 


9.8 - 11.7 


70 



a Samples of star-forming galaxies studied in this paper. 

lection crite ria are based on the observed U n GR colors 
of galaxies. IReddv et al.l (|2006l ) provided spectroscopic 
observations for candidates b righter than Rab = 25.5 , 
photometrically selected from ISteidel et all ((20031 12004) . 
For this paper, 41 BM/BX galaxies to Kab — 23 are 
selected from their sample with spectroscopic redshift 
z > 1.4. 

Galaxies analogous to the LBGs at 0.6 < z < 
1.4 (hereafter GALEX/LBGs) are selected with the 
GALEX/HST (FUV - NUV) AB versus (NUV - 
F435W)ab color-co lor diagram. Follow ing the selec- 
tion criteria used in iBarger et al.l (|2008 ii) and requiring 
z > 0.6, we include 105 GALEX/LBGs with K AB < 23. 

In addition to the UV-bright galaxies in GOODS-N 
field, we also study sizes of the Sub-millimetre galax- 
ies (SMGs) as a population of star forming galaxies at 
high redshifts. The subsample of SMG s fo r our stud 
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signed to find actively star forming galaxies at redshifts 
1.4 < z < 2.5 with similar SEDs to the LBGs. The se- 



ies is d rawn fromlChapma net al.l (|2005| ) and 
(2006). The catalog of SMGs provided by 
(2006) contains 35 candidates from GOODS-N field with 
21 secure optical counterparts. Of these, only 8 galaxies 
in our observed images have spec troscopic redshif t s. We 
further include 6 SMGs from the lChapman et al.l (|2005[ ) 
HDF-North study. Therefore, our final SMG sample 
studied here contains 14 sources spanning 0.5 < z < 3. 

Besides these two populations, star forming galaxies 
at z > 1.4 can als o be identified using BzK technique 
(jDaddi et al.l [2004). This technique is designed to cull 
galaxies from i^-selected samples. Star forming BzK 
galaxies (sBzK) tend to be more massive and have higher 
reddening than UV-selected galaxies. With this method, 
70 sBzK galaxies at z > 1.4 and Kab < 23 are selected 
using the criteria (z — K) — (B — z) > —0.2. We further 
use this sample to compare properties of star forming 
galaxy populations at high redshifts. The sBzK selec- 
tion criteria can identifies ~ 90% of our BM/BX galax- 
ies. The large overlap is not surprising since both our 
UV-selected and sBzK samples consist of mostly mas- 
sive galaxies (due to the K magnitude limit). Table Q] 
lists samples of star forming galaxies used in this paper. 

6. RESULTS 

6.1. Size Evolution 

In Figure [3] we show the size distributions for UV- 
bright galaxies between 10 < log(M^/M Q ) < 11. The 
blue histogram shows the measured half-light radii of 
GALEX/LBGs at z ~ 1 while the distribution of sizes 
measured for BM/BX galaxies (z ~ 2) are shown as a 
green histogram. It can be seen from this plot that UV- 
bright galaxies at z ~ 2 are smaller compared to similar 
galaxies at z ~ 1. Specifically, the median half light 
radii of BM/BX galaxies in this mass range is 2.68 ±0.19 
kpc significantly smaller than the median effective radii 
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Fig. 3. — Size distribution of UV-bright galaxies with stellar 
mass 10 < log(M-*/Mo) < 11 at z ~ 1 (blue histogram) and 
z ~ 2 (green histogram). The half light radii of UV-bright galax- 
ies (GALEX(LBGs) ) at z ~ 1 are larger than BM/BX galaxies at 



of GALEX/LBGs (4.42 ± 0.52 kpc). This means that 
sizes of UV-bright galaxies evolve by a median factor of 
0.60 ± 0.08 between z ~ 2 to z ~ 1. It is harder to mea- 
sure the evolution for z ~ 3 to z ~ 2, as we have 4 LBGs 
within mass range of 10 < log(M ie / Mq) < 11 at z ~ 3 
in our sample. The median half-light radius of these 4 
massive galaxies is 2.22 ± 0.61 kpc, 0.82 ± 0.23 smaller 
than z ~ 2 BM/BX galaxies. We perform a power-law 
fit on the size evolution and find r e oc (1 + z ) _111±013 
over the range 0.6 < z < 3.5. 

The sizc-rcdshift relation for star-forming galaxies with 
spectroscopic redshifts in the GOODS North field is 
shown in Figure 01 Galaxies are split into the four differ- 
ent stellar mass bins shown in the figure. In each panel, 
the gray small symbols are normal star forming galax- 
ies with spectroscopic redshifts selected to have specific 
star- formation rates log sSFR > —10. Color symbols 
represent our samples of UV-bright galaxies and sub- 
millimetre galaxies (SMGs). LBGs (2.7 < z < 3.5) 
are shown as dark-green triangles and BM/BX galax- 
ies (1.4 < z < 2.7) are plotted as green squares. 
GALEX/LBGs at relatively lower redshifts (0.6 < z < 
1.4) are plotted in different panels as blue stars. The 
black sym b ols at z ~ are star forming galaxies from 
IGuo et all ((2001 from SDSS. 

As Figure @] shows, sizes of UV-bright galaxies at a 
fixed stellar mass increase towards lower redshifts. There 
is also a trend that the size evolution maybe stronger 
for the most massive galaxies. For example, UV-bright 
galaxies with stellar masses 9 < Zog(M»/M Q ) < 10 
evolve by a median factor of approximately 1.15 ± 0.25 
from z ~ 2 to z ~ 0.8. However, galaxies with masses 
10 < Zog(M*/M Q ) < 10.5 grow by a median factor 
of approximately 1.82 ± 0.25. At higher stellar mass 
bin (10.5 < log(M*J Mq) < 11) this evolution is even 
stronger: a factor of 2.05 ± 1.26 over the same redshift 
range. However, selection effects may influence this trend 
(especially in the low-mass bin) and more complete sam- 
ples are needed. 

SMGs are also shown as pink diamonds in Figure |4j 



they are discussed further in section 6.4. Galaxies which 
have X-ray luminosities > 10 42 ergs s _1 (i.e. potential 
AGN hosts) are marked with open circles. These galaxies 
have sizes similar to the others. Since AGN is point 
like, one would expect galaxies with AGN emission to 
have smaller half-light radii. Their "normal" rest frame 
optical sizes are thus perhaps not strongly influenced. 
Nonetheless, due to possible effects on the size and mass 
estimates, X-ray detected galaxies should be considered 
low-confidence points. 

6.2. Stellar Mass-Size Relation 

We now investigate how mass-size distribution of the 
UV-bright galaxies compares to other galaxy popula- 
tions. We further check if the relation between stellar 
mass and size of these galaxies exists to high redshift. 

The stellar mass-size distributions of our samples are 
shown in FigureO We have divided our sample into three 
redshift bins, 0.5 < z < 1.5, 1.5 < z < 2.5, 2.5 < z < 3.5. 
UV-bright galaxies are color coded the same as Figure |4j 
Solid and do tted lines in each panel show the mass-size 
relation from lShen et al.1 (|2003fl for star forming galaxies 
and early-type galaxies at z ~ 0, respectively. The gray 
symbols are the if -selected sample o f galaxies in CDF - 
South with stellar masses > 1O 9 8 M (jFranx et al. 2008). 
Comparing our results with the galaxies from CDF-S al- 
lows us to see where our galaxies lie relative to a purely 
mass-selected sample. Brown solid lines in the middle 
of each panel show the median effective radii in narrow 
mass bins for all galaxies in CDF-S. As can be seen, at 
z ~ 1 and z ~ 2, UV-bright galaxies have larger effective 
radii compared to all galaxies at a given stellar mass. 
For galaxies between 10 10 -10 1:l Mq, the median sizes dif- 
fer by a factor of 1.89 ± 0.23 and 1.12 ± 0.09 at z ~ 1 
and 2 respectively. As the spectroscopic sample contains 
only 4 UV-bright galaxies in this mass range at z <~ 3, 
the difference with CDF-S is not well constrained at this 
redshift (1.27 ±0.36). 

We further investigate how the size mass distribution 
of GOODS-N UV-bright galaxies compares to that of 
UV-bright galaxies in CDF-S. This would help to verify 
if the previously rep orted size-mass r elation for a large 
photometric sample (jFranx et al.ll2008T ) is consistent with 
our spectroscopic sample. We define analogues of the 
z ~ 1,2 and 3 LBGs in CDF-S by their optical colors 
and magnitudes. Specifically, we selected the analogues 
of the z ~ 1 GALEX/LBGs by requiring (B — V) < 0.6 
and R < 23.8. This selection is based on the locus of 
the GOODS-N GALEX/LBGs in this color-magnitude 
plane. The analogues of BM/BX galaxies in CDF-S are 
similarly selected with R < 25.1 and [B — V) < 0.8 and 
Zphot = 1.5 — 2.5. We note that the R band magnitude in 
both surveys is estimated by the average values of V and 
/ band magnitudes. We verified that the spectroscopic 
samples are effectively selected by these criteria. In order 
to identify analogues UV-bright objects in CDF-S at z ~ 
2.5 - 3.5, we find the location of LBGs in GOODS-N 
on the color-color (B — V) versus (V — I) diagram and 
apply the following criteria to the CDF-S galaxies; (1) 
0.5 < (B — V) < 1.2 for (V - I) < 0.35 ;(2) (B-V)> 
3 X (V -I) + 0.5 and (B — V)< 1.2 for 0.35 < (V - 1) < 
0.5. We note that these selection criteria do not impose a 
"U-drop out" criteria, but do select galaxies with similar 
brightness and UV slope as the usual LBG methods. 
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Fig. 4. — Evolution of sizes of galaxies as a function of redshift. Each panel shows the size evolution for a narrow stellar mass bin. The 
gray symbols are star forming galaxies with spectroscopic redshifts in GOODS-N field. Blue stars are GALEX/LBGs (z ~ 1) while the 
green squares and green triangles represent BM/BX galaxies (z ~ 2) and LBGs (z ~ 3)respectively. Submm galaxies are also shown as 
pink diamonds. Galaxies with X-ray l umin osities more than 10 42 ergs s~ 1 are marked by open circles. Black symbols at z ~ are SDSS 
star forming galaxies from Guo et al. (200(J). As this plot illustrates, the half light radii of UV-bright galaxies evolve with redshift and this 
evolution is faster for high mass galaxies. 



The comparison of size-mass distributions for UV- 
bright galaxies from both fields is shown in Figure [5] In 
each panel, the color symbols are the median log effective 
radii of UV-bright galaxies in narrow mass bins, with the 
blue and green symbols representing UV-bright galaxies 
in GOODS-N and the brown symbols are the ones in 
CDF-S. Error bars show the one a dispersion. UV-bright 
samples from both fields follow consistent size-mass re- 
lation, further confirming that UV-bright galaxies are 
on average larger than overall galaxies at a fixed stellar 
mass. We note that at high redshift (z ~ 2.5 — 3.5) the 



comparison is weak because of the limited number LBGs 
in GOODS-N. 

The UV-bright galaxies show a weak stellar mass size 
relation at z = 1 and z — 2 in Figure [5j To quan- 
tify this relation, we fitted a power law function of the 
form r e oc M a to the individual UV-bright galaxies in 
each redshift bin. The fitting results are plotted as color 
dashed-dotted lines in Figure [51 At z ~ 1, the size of 
GALEX/LBGs in GOODS-N scales with stellar mass 
as r e oc M 019±0 05 and the BM/BX galaxies at z - 2 
have r e oc jV^°-30±o.06^ uncer ^ am tj es W ere estimated 
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Fig. 5. — Stella r mass size distribution for UV-brigbt galaxies in different redshift bins compared to the galaxies from CDF-south 
IjFranx ct al. 2008]) (gray dots). The color symbols represented here are the same as Fig. [4] The sol id and dotted black lines are the 
size-mass relations for star forming and quiescent galaxies respectively at z ~ from Shen ct al. (2003). The solid brown lines show the 
median sizes in narrow mass bins for all galaxies from CDF-S. The UV bright galaxies are in general larger than normal field galaxies at 
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Best fits power law parameter for the stellar mass 
relation 



Fig. 6. — Comparison of the stellar mass-size relation for the UV-bright galaxies from GOO DS-N field a n d CD F-South in three redshift 
bins. The black solid lines are the size-mass relation for star forming galaxies at z ~ from Shen et al. (2003). The color symbols are 
median of UV-bright galaxies in narrow mass bins and the error bars show one a (68%) dispersion. The dashed-dotted lines are the best 
power law fits using r e oc M a to the individual UV-bright galaxies. The dotted lines are the median sizes in narrow mass bins for all 
galaxies from CDF-S (Same as Fig [5}. This plot shows that the derived half-light radii of UV-bright galaxies from both fields are in place 
and agreement with each other up to z ~ 3. 

using bootstrap resampling. The results are compara- 
ble to the r elatio n for late-type galaxies at z ~ from 
iShen et all (|2003f ). They found r e oc M 015 for low-mass 
(log(M) < 10.6) late-type galaxies and steeper relation 
for high mass late type galaxies (r e oc M 0A ). In the high- 
est redshift bin (z ~ 3), we use UV-bright galaxies from 
CDF-S to find the best fit (a = 0.32 ±0.06). Tabledlists 
the best-fit power law parameter to the mass size relation 
of UV-bright galaxies. The slopes of the mass-size rela- 
tion at different redshift bins are consistent and there is 
no significant evolution. Our results confirm the persis- 
tence of the size-mass relation for star forming galaxies 
up to high redshift. 

It is worth checking whether or not other star form- 
ing populations at high redshift have the same size-mass 
relation. Therefore, we compare star forming BzK galax- 
ies (sBzK) with BM/BX galaxies. The size distribu- 
tion as a function of stellar mass is shown in Figure [7] 



Sample 


Redshift 




a 




GALEX/LBG 


0.6 < 


z 


< 1.4 


0.19 


± 


0.05 


BM/BX 


1.4 < 


z 


< 2.7 


0.30 


± 


0.06 


CDFs-UV Bright 


2.5 < 


z 


< 3.5 


0.32 


± 


0.06 


sBzK 


1.4 < 


z 


< 3.2 


0.31 


± 


0.09 



a Note. Power law parameter a is defined as r e oc M a . 

with the sBzK galaxies plotted as gray filled circles. The 
median of log effective radii of sBzK galaxies are over- 
plotted as red squares. BM/BX galaxies and SMGs are 
marked with the green and pink circles respectively and 
the green dashed-dotted line shows the mass-size rela- 
tion for BM/BX galaxies. As one can see, the mass- 
size relation for sBzK galaxies is comparable to the one 
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Fig. 7. — Stellar mass size relation for star forming sBzK galaxies 
(gray circles). The BM/BX galaxies and SMGs are marked with 
green and pink circles respectively. The red squares are the median 
for the mass-size relation of sBzK sample with one a dispersion. 
The green dashed-dotted line is the best fit to the BM/BX galax- 
ies. The sBzK galaxies have a similar size-mass relation to the 
BM/BX galaxies. The solid line s hows the mass-size relation for 
star-forming galaxies at z ~ from Shcn ct al. (2003). 

for BM/BX galaxies. The average offset of median half 
light radii of sBzK galaxies to the BM/BX mass-size re- 
lation is 0.037 dex. The effective radii of sBzK galaxies 
scale with stellar mass as r e oc M - 31±om close to the 
size-mass relation for UV-bright galaxies. This similarity 
could be due to the significant overlap bet ween the UV- 
select ed galaxies and sBzK galaxies (e.g. iReddv et alJ 
12005( 1 . We note that stellar masses can have significant 
uncertainties depending on the methods used to calculate 
the m. We verify that by using the stellar mass estimates 
from lReddv et al. (2006}) for UV-bright galaxies, the size- 
mass distribution of these galaxies covers the same region 
in the size-mass plane (see e.g., middle panel of Figure 
[5|). Hence, this shows that our results are robust against 
using different estimates of stellar masses. 

6.3. Color - Surface Density 

One of the surprising results of lFranx et al.l (|2008f l was 
the tight correlation between stellar mass surface density 
(M*/fg) and color of galaxies to z ~ 3. They showed 
that bluer galaxies have lower surface densities than red 
quiescent ones. They indicated that the color of galax- 
ies correlates more fundamentally with the stellar mass 
surface density than the mass. 

We show the tight relation between color and stellar 
mass surface density for both star-forming and quiescent 
galaxies with M & ~ 10 10 - 10 11 in Figure M In the left 
panel the observed color (Veoe — zsso) is plotted versus the 
stellar mass surface density for galaxies at redshift ~ f . 
The black and gray symbols are galaxies from GOODS-N 
and CDF-S fields, respectively. The blue stars and brown 
circles indicate the GALEX (LBGs) and their analogues 
from CDFS respectively. In the right panel, the rela- 
tion between observed (zsso — K) and surface density 
for galaxies at higher redshift ~ 2 is illustrated , with 



BM/BX galaxies in GOODS-N shown by green squares 
and their analogues in CDF-S as brown circles. As can 
be seen in this Figure, the blue star forming galaxies have 
lower surface density than red galaxies at both z ~ 1 and 
z ~ 2, though the correlation between color and surface 
density is more tight at z ~ 1 than z ~ 2. In general, 
using sample of galaxies with spectroscopic redshifts con- 
firms that the color of galaxies are tightly correlates with 
surface density and this relation holds out to high red- 
shifts. 

6.4. Sizes of Submm Galaxies 

Since GOODS-N includes 14 SMGs with redshifts, we 
compare their mass and sizes to the UV-selected galaxies 
at z ~ 2. As Figure |4] illustrates, the optical rest-frame 
sizes of the SMGs (pink diamonds) and UV-selected 
galaxies of similar mass are comparable. The median half 
light radii of all SMGs is 2.90 ±0.45 kpc which is in good 
agreement with the median optical rest-frame sizes of the 
BM/BX galaxies over the whole mass range (2.68 ± 0.25 
kpc). The median effective radii of SMGs with stellar 
masses of 10 10 — 10 11 is 2.65 ± 0.56 which is similar to 
the BM/BX of at the same mass range (2.68±0.19). The 
errors are computed using bootstrap resampling. 

We further compare the size distributions of sub-mm 
galaxies and BM/BX galaxies in Figure H] In the 
top panel, the histograms show the size distribution of 
BM/BX galaxies and SMGs. In the bottom panel, the 
normalized cumulative distribution functions are shown 
for both samples. As this plot shows, the size distribu- 
tions of SMGs are comparable with z ~ 2 UV-bright 
galaxies and the SMGs follow the same size distribu- 
tion as BM/BX galaxies in the optical rest- frame with 
SMGs being slightly larger. According to a KS test, 
the probability that two distributions being the same is 
72 percent, suggesting no significant difference between 
the rest-frame optical sizes of the SMGs and UV-bright 
galaxies. 

7. DISCUSSION 
7.1. The growth of star-forming galaxies 

We study the size evolution of galaxies with spectro- 
scopic redshifts between z ~ 0.5 — 3.5. Our results are 
summarized in Figure QUI As this plot illustrates, at 
each epoch, the median sizes of UV-bright galaxies (color 
squares) with stellar masses 10 < log(M*/ M©) < 11 are 
larger (by a median factor of 0.45 ± 0.09) than quiescent 
galaxies in a similar mass range selected from CDF-S (red 
circles). The red tri angle is the median s i zes of q uiescent 
galaxies studied in Ivan Dokkum et alJ (|2008b( l (VD08 
hearafter) with a median stellar mass of 1.7 X 10 n M Q 
and median redshift of 2.3. Table |3] lists the median half 
light radii for UV-bright and quiescent galaxies seen in 
Figure[TUl The growth of UV-bright galaxies with time is 
also illustrated in this figure. The dashed line shows that 
UV-bright galaxies scaled up their sizes towards lower 
redshifts as (1 + z) - 1 - 1 ^ 013 . 

Pre vious studies (jFranx et al.l 120081 : IKauffmann et al.1 
12003( 1 reported that there is a correlation between color 
and stellar mass surface density of galaxies at both low 
and high redshifts. However, the results based on photo- 
metric redshifts can be uncertain especially for star form- 
ing galaxies. By means of our large sample of galaxies 
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Fig. 9. — Top panel: Histograms show the distribution of sizes 
of Submm galaxies (solid line) and BM/BX galaxies (dotted lines). 
The gray region indicates apparent effective radii below 0.1" (~ 1 
kpc at 2 ~ 2) which size measurements have large uncertainties. 
Lower panel: Cumulative distribution function for SMGs (pink 
dimonds) and BM/BX galaxies (green squares). As the plot shows, 
the distributions of half-light radii of these two populations are 
comparable with SMGs being slightly larger. 

with secure redshifts, we have confirmed that the tight 
correlation between color and stellar mass surface den- 
sity of galaxies that holds out to high redshifts; as the 
stellar mass densities of blue star forming galaxies are 



J Sample 


Redshift 


Size 


median(log(M* / Mq) 


GALEX/LBG 


0.6 < z < 1.4 


4.42 ±0.52 


10.2 ± 0.3 


; BM/BX 


1.4 < z < 2.7 


2.68 ±0.19 


10.4 ±0.2 


: LBG 


2.7 < z < 3.5 


2.22 ±0.61 


10.6 ± 0.2 


: CDFS 


0.5 < z < 1.5 


2.33 ± 0.1 


10.4 ±0.3 


~ CDFS 


1.5 < z < 2.5 


2.40 ±0.13 


10.4 ±0.3 


: CDFS 


2.5 < z < 3.5 


1.73 ±0.18 


10.4 ±0.3 


Quiescent(CDFS) 


0.5 < z < 1.5 


1.32 ±0.07 


10.5 ± 0.3 


; Quiescent(CDFS) 


1.5 < z < 2 


1.12 ±0.32 


10.7 ±0.2 


- Quiescent(VD08) 


2 < z < 2.5 


0.9 ±0.21 


11.23 


: SMG 


0.5 < z < 3. 


2.90 ±0.45 


10.8 ± 0.5 



a Notes. The median sizes of samples are for galaxies with stellar 
masses between 10 10 — IO ^Mq except for the quiescent(VD08) 
Ivan Dokkum et al.ll2008"EI) and SMG galaxies. The CDFS sample 
is all galaxies at the same stellar mass range in CDF-S field. 

smaller than those of red quiescent ones. This verifies 
that the galaxies with higher specific star formation rate 
are larger than the ones with lower specific star formation 
rates. 

We have also explored the stellar mass-size distribution 
for galaxies in our sample. We have confirmed that there 
is a relation between stellar mass and size for UV-bright 
galaxies at 0.5 < z < 3.5. The relations are consistent 
for both spectroscopic and photometric redshift samples. 
Our results verify that there is no significant evolution 
of the slope of the size-mass relation to redshift ~ 3.5. 
The evolution of the relation from z = 1 to z = is 
not well established, as it is not straight forward to se- 
lect Lyman Br eak galaxies at z ~ 0. Some authors find 
evolution (e.g. iFranx et al.l 120081: iWilliams etldl 12010ft 
in the stellar mass-size relation, however. iBarden et al.l 
(2005) have reported that the stellar mass-size relation 
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Fig. 10. — The median sizes of UV-bright galaxies (squares) 
as a function of redshifts for galaxies with stellar masses 10 10 < 
M/M© < 10 11 in GOODS-N field. The red filled circles are qui- 
escent galaxies from CDF-S study w ith similar mass range and 
red triangle is quiescent sample from van Dokkum ct al. (2008b) 
with median stellar masses of 1.7 X lO^M©. The dashed line 
shows the best fitting size evolution to the UV-bright galaxies 
(r e oc (1 + 2)-i ii±o i»). The plot indicates that the UV-selected 
galaxies are larger than quiescent galaxies and their sizes evolve 
with redshift. 

for disk galaxies remains constant from z <~ 1 to the 
present. Selection effects likely play a role. Comparing 
to the theoretica l work , the size evolution predicted by 
iSomerville et al.l (pO OS) is somewhat stronger than the 
evolution we measured here. 

Star forming galaxies at high redshifts can be identi- 
fied by means of different methods of sample selection 
and can be assigned as different populations of galaxies 
(e.g, sub-mm or sBzK) and therefore studied separately. 
However, many of these galaxy populations may over- 
lap, hence it worth checking out if they have comparable 
properties. We therefore, compare sizes of UV-bright 
galaxies to other types of star-forming galaxy popula- 
tions at high redshift, i.e. sBzK galaxies and SMGs. We 
show that these two populations have comparable sizes 
to the UV-bright galaxies. This suggests that in gen- 
eral, star forming galaxies have larger sizes than quies- 
cent ones without regard to their type. However, sample 
selection might alter the results. We return to this pos- 
sibility later, in the section 7.3. 

7.2. Comparison to previous work 

The strong growth of both star forming and qui- 
escent galaxies has been shown by s everal authors 
dTruii lo et al.1 l2006ri IF ranx et al.l [200l iWilliams et all 
l2010t iToft et al.l 12007ft . The size growth of star form- 
ing galaxies studied here going as (1 + z) a where a = 
— 1.11 ± 0.13, agrees well compared to previous stud- 
ies where the size evolve s with a ~ — 1 . For exam- 
ple, iBouwens et ail (p006l ) found a = -1.1 ± 0.3 for a 
sample of luminosity-selected dropout galax i es be tween 
z ~ 6 and z ~ 2.5 while iDahlen et all (120071 ) find 
a = —1.10 ±0.07 for lumin osity-selected d i sk gal axies 
over the range z ~ 1.1 - 2.2. IWilliams et al.l ()2010D also 
find a similar slope for star forming galaxies in UDS at 



z ~ 0.5- 2. 

The similarity between sizes of UV-bright galaxies and 
star forming BzK galaxies is also show n in our study. 
The tw o populations are also compared bv lOverzier et all 
(2010). They showed that rest-frame optical sizes of 
sBzK galaxies are somewhat larger than UV-selected 
galaxies. Since BzK galaxies are K-selected and thus 
tend to represent massive galaxies, it is not surprising 
that they appear larger. Hence, the sizes must be cor- 
rected for differing stellar masses for an accurate com- 
parison to be made. 

Our results also show that the SMGs have a median 
half light radii of 2.90 ± 0.45 kpc comparable to the rest- 
frame optical sizes of the BM/BX gala xies. The result 
is in a greement with a recent study by [Swinban k et al.l 
(2010) where they use deep HST I and if-band imag- 
ing and show that the rest-frame optical sizes of the 
SMGs and UV-bright galaxies in GOODS-N are compa- 
rable. The typical half-light radii of the SMGs in the H- 
band (rest-frame optical) measured by ISwinban k et al.l 
(2010) is 2 .8 ± 0.4 which is co nsistent with our mea- 
surements. lAlmaini et al.l (|2005l ) also compared sizes of 
submm galaxies with Lyman-break galaxies at rest-frame 
optical and found no clear difference between sizes of the 
SMGs and LBGs. It is worth noting that SMGs have 
larger median stellar masses compared to BM /BX galax- 
ies (also mentioned in ISwinbank et ail (|2010T )). This sug- 
gests higher stellar mass densities for SMGs compared to 
UV bright galaxies. However, the similarity in sizes be- 
tween these two galaxy populations at high z does not 
by itself allow conclusive connections to be drawn be- 
tween the two populations. The SMGs that are faint or 
lacking emission lines could be missing from our spectro- 
scopic sample and they may have dramatically different 
mass or size properties. Therefore, selection effect could 
bias these results. Spectroscopic observati on in near- 
IR mig ht help to solve this problem; e.g., iKriek et al.l 
(2009a) find a median size of 2.8 kpc for a small sam- 
ple of near-IR spectroscopically confirmed star forming 
galaxies at z ~ 2.3. Indeed, our measured size evolu- 
tion is consistent with their spectroscopic sample; how- 
ever, near-IR spectroscopy over a large redshift range is 
needed to definitively quantify the size evolution. 

7.3. caveats 

With this large spectroscopic sample of high redshift 
galaxies, we have removed an important source of 
uncertainty in the stellar mass size relation and verified 
the results from larger samples based on photometric 
redshifts. However, the spectroscopic catalog we used 
in this study consists of many different spectroscopic 
surveys which made it inhomogeneous. Galaxies in 
this catalog are mostly relatively unobscured star 
forming galaxies with emission lines, comparing these 
to quiescent galaxies can still be a problem because 
spectroscopy of high redshift red galaxies is difficult. 
Therefore, using a larger and homogeneous sample of 
high redshift galaxies with secure redshift will allow 
greatly improved constraints on their evolution. New 
ground and space-based IR multi-object spectrographs 
(WFC3 grism, MMIRS, MOSFIRE) will allow inves- 
tigations of more complete samples of galaxies and 
studies the correlation between galaxy properties at 
higher redshifts. Dust-obscured starbursts are also 
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likely to be missed in our UV-bright sample, hence NIR 
spectroscopy will also permit studies of this important 
population. 



8. SUMMARY 

We have performed the first size evolution study of 
UV-bright galaxies with secure redshifts in the GOODS- 
North field. We derive half light radii of galaxies at a 
wide range of redshifts (up to z ~ 3.5). Using sample of 
galaxies with secure redshifts, we have quantified the size 
evolution of galaxies without the potential uncertainty of 
photometric redshifts and confirm the previous studies 
that galaxies of similar mass were generally smaller at 
past and their sizes evolve with redshift. Specifically, we 
find that: 

• UV-bright galaxies evolve strongly as (1 + 
^-i.ii±o.i3 connrrmn g prev ious studies based on 
photo metric redshifts (e.g. iTruiillo et ail [2 006b. 
120071 : iFranx et al.ll2008t iWilliams et al.ll2010f ). 

• UV-bright galaxies are significantly larger than qui- 



escent galaxies. At the same mass, the median dif- 
ference is 0.45 ± 0.09 dex. 

• The LBG, BM/BX and GALEX/LBG samples 
show smooth evolution with redshift, indicating 
that these techniques indeed select similar UV- 
bright galaxies at different redshifts. 

• The SMGs have half-light radii similar to UV- 
bright galaxies of the same mass. 
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